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' Abstract 

■ The partial decay rate T(Z — ► 66) is significantly influenced by the 
mass of the top quark due to electroweak radiative corrections. The 

■ leading ~ mf and the next-to-leading contribution ~ In mf are known 
"c^ to be numerically of similar size. In this work we calculate the QCD 

corrections to the logarithmic correction using the heavy top mass 
expansion. The 0(aa s lnm^) corrections are of the same order as the 
^ ■ QCD corrections to the quadratic top mass term, but of different sign. 
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1 Introduction 



Although the direct observation of the top quark is out of range for the experi- 
ments at LEP, several observables are affected by the top quark through virtual 
states in higher order radiative corrections. 

High precision measurements and the comparison of these quantities with 
the theoretical predictions allow to extract bounds on the top mass. Present 
analysis of e + e~ collisions at the Z peak estimates the top mass in the range m t = 
173+]Jt2i GeV [|IJ, which is in agreement with top masses of m t = 174 ± lO^ 
GeV from pp collisions at TEVATRON 0. Of particular interest for deducing 
the limits on m t from LEP data is the partial decay rate YiZ — > bb). On the one 
hand this quantity exhibits a strong sensitivity on the top mass @ as the leading 
term of the electroweak corrections is quadratic in m t and the next-to-leading 
logarithmic contribution In m 2 is numerically of the same size. On the other hand 
the already small uncertainty of the measurement of R h i = YiZ — > bb)/T(Z — > 
hadrons) = 0.2192 ± 0.0018 U is expected to be reduced below 1% in the future. 
As a consequence the determination of QCD corrections to the electroweak one 
loop result became increasingly important. 

For the quadratic top mass contribution these 0(aa s m 2 ) corrections were 
calculated by four independent groups [|, |§, §|, [3 • in this work we present the QCD 
corrections of order 0(aa s him 2 ) to the next-to-leading top mass contribution. 

Our calculation is performed using dimensional regularisation in the MS- 
scheme with anticommuting 75 in the 't Hooft-Feynman gauge. As in our previous 
work |7j we employ the hard mass procedure || ||, [U], [UJ to derive the next-to- 
leading order in the inverse top mass expansion. In order to handle large expres- 
sions during this expansion we use the algebraic manipulation program FORM 
[fT^I - Massless multiloop integrals are evaluated with the help of the software 
package MINCER |[J. 

If we consider the properly resummed propagator of the Z boson we get in 
the case q 2 m M§ for the transverse part 

n^fn 2 ) = -- l - ^ m 

i 1 + U'(Ml) M 2 n 2 »Im n(Mf )M Z W 

v zj m z -q T+n 7 (Mj) 

where Mz is the renormalized Z boson mass and II (q 2 ) is the transverse part 
of the polarisation tensor H^q 2 ). One can get II(g 2 ) through contracting with 
l/(D — l)(g^ v — l^ilu/l 2 ) where D = 4 — 2e is the space-time dimension. n'(M|) 
is defined as d Rell(q 2 )/dq 2 also evaluated for q 2 = Mf. Thus we can write the 
partial decay rate YiZ — > bb) in the following form 



T(Z -> bb) 



T (v 2 + a 2 ) + Ar 



1 + ir(Mf) 
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with v = — 1 + fsvp,o = — 1 an d To = l4$>Cw s w- ^ contains all vertex 
corrections and the 7Z-mixing. Expressing the result in terms of Gp (instead of 
a) introduces a factor (1 — Ar). The above equation can be written in the form 

■GpMl 



T(Z -> bb) 



" " z (v 2 + a 2 )+AT 



1 - Ar 

i + n'(M|; 



(3) 



The top mass dependence of the last term is of universal nature. This contribution 
- as well as the universal part of Ar — can be expressed through the universal 
p-parameter |T3, 15]. In this work we calculate the non- universal part of order 
0(aa s him 2 ). 



The outline of the paper is as follows: Section 2 contains further details of the 
calculation and in section 3 we present the results in the MS and OS-scheme. 



2 Calculation of the Order 0(aa s \nm^ 

In the decay rate F(Z — > bb) the top mass dependence is due to the appearance 
of the top quark as a virtual particle, i.e. it is induced by t — 6-transitions due 
to the exchange of charged Higgs ghosts $ =t or W bosons. One can distinguish 
between seven classes of diagrams as listed in figure [I| from which the imaginary 
part has to be calculated. QCD corrections are obtained by attaching gluon lines 
in all possible ways. This results in 58 topologically different diagrams. 

We perform our calculation in the large top mass limit m t — > oo, i.e. we 
apply an asymptotic expansion in the inverse heavy mass l/m t || |§ [H| [TT 



In practice one isolates all possible hard subgraphs and expands them w.r.t. the 
small masses and (external) momenta. Afterwards this expansion is inserted as 
an effective vertex into the remaining diagram. 

The hard mass procedure was already used in our previous work |7| to evaluate 
the 0(aa s m 2 ) corrections. In this order only diagrams with Higgs ghost exchange 



Figure 1: The seven distinct classes of electroweak diagrams, which contain the 
top quark. Dashed lines: Higgs ghost; wavy (internal) lines: W boson; thin lines: 
b quark; thick lines: t quark. 
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contributed. The calculation of the next-to-leading 0(aa s himf) corrections is 
extended in three different ways. First, the diagrams with W boson exchange 
need to be taken into account. Second, additional hard subgraphs of the Higgs 
ghost diagrams contribute to the considered order. Third, the subgraphs already 
considered in |7j] must be expanded to the next higher order. 

We find that the next-to-leading order results are finite separately for each 
class of diagrams which contain a Higgs ghost. Third order poles 1/e 3 compensate 
in the sum of all possible hard subgraphs of a diagram and second order poles 
1/e 2 cancel after replacing the bare top mass by the MS renormalized one 

The remaining simple divergence drops out in the imaginary part of the diagram. 

The class of W diagrams and their renormalized contribution RU^ w \q 2 ) re- 
main to be calculated. Here R denotes the so-called Bogoliubov-Parasiuk R- 
operation [|1(J adopted to the minimal subtraction scheme [jT^] . This procedure 
determines the finite renormalized value of a given regularized Feynman integral. 
It works on a graph-by-graph basis and removes all sub divergencies together with 
the overall UV divergence in a way compatible with adding local counterterms to 
the Lagrangian. 

We apply this prescription to subtract all sub divergencies of each W diagram, 
whereas the overall divergence again is eliminated by taking the imaginary part. 

The pole part A (7) of a divergent subgraph 7 of a Feynman diagram F is 
defined as the overall divergence of the subgraph 7 and can be written as a 
polynomial in possible masses and external momenta of 7 with coefficients being 
pure poles in e. Since the e poles themselves are independent of masses and 
momenta, their determination can be simplified by setting masses or momenta 
zero. However, care must be taken to guaranty that no spurious IR singularities 
are introduced this way. The subtracted Feynman integral A (7)^/7) is obtained 
by replacing the subgraph 7 through A (7). For the calculation of the remaining 
integrals (r/7) the heavy mass expansion can be employed. 

There are several checks for the correctness of our result. The first one is the 
fact that after renomalization the seven classes of diagrams are separately finite. 

The second check is the gauge invariance with respect to QCD. The calculation 
is done with an arbitrary QCD gauge parameter £s which is introduced through 
the gluon propagator. As expected each class is separately independent of £5 if 
the sum of all diagrams is taken. 

The third check is the invariance of the gauge parameter £iy from the elek- 
troweak theory. £vk appears through the propagator of the W boson and the $ ± . 
The total sum of all diagrams is indeed independent of 
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3 Results and Discussion 



In order to calculate explicitely the 0(aa s him 2 ) correction to the Z —>■ bb vertex, 
we have employed the hard mass procedure as an expansion up to the next-to- 
leading order in the 1/rrit series. As a consequence the leading order result is 
automatically reproduced. Furthermore it became necessary to repeat the calcu- 
lation up to next-to-leading order for the case without QCD corrections (see fll8| ) 
since they induce corrections in first order a s through top mass renormalization. 

Let us first recall this purely electroweak result for the non-universal decay 
rate. 




The logarithms have their origins in different classes of diagrams. From the 
diagrams with Higgs ghost exchange logarithms In [i 2 /fn\ and ln// 2 /Mjy occur 
such that in the combination the \i drops out and In M^/m 2 are left. The In /x 2 /m 2 
only comes from the W graphs. Choosing the scale fi 2 = reproduces the 
result known from ||. Since we are only interested in top mass effects, i.e. in 
the difference of a heavy and a light top quark, we are allowed to set [x 2 = 
because the difference T(Z — > bb) hcavytop — T(Z — > bb) Ughttop is independent of fi 
and only contains logarithms of the form In m 2 /M^ v . 

Including first order QCD corrections we obtain the non-universal part of Ar 
in the MS scheme: 



G F _ 2 2 ^ 
V2n 2 V 3 



7T 



MS 
rem 



+ Ar" 
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Ar rom contains the constant terms of the non-universal diagrams containing 
the top quark and also terms from graphs without top quark which we have not 
calculated. Via the relation 



m(fi 2 



m Q s 



1 _ «s 

7T 



4 , 
- + ln 



m 



OS 



we transform this result into the OS-scheme: 



Ar 



non— univ. 



Gf I ' 2 2 \ a s (jj 



*-* °V2^\ 3 W J 7T 

f u 2 

Mfj In ^ (-3 + 6^-35^) 



7 , Mf / . . 2 
T2" 



+ Q-t* + 34,1+^1- 



rem 



■«w G F / 2 2 \ a. 



mi V 3 



em 



As a consequence of dimensional regularisation the renormalization scale // 
appears in the calculation. For the diagrams with Higgs ghost exchange again the 
logarithms In ^ 2 /m 2 , In /i 2 /M^ and ln/i 2 /M§ are combined such that /i 2 disap- 
pears. Since these logarithms are connected with the pole structure of the result, 
this cancellation is expected from the fmiteness of these classes of diagrams. 

For the classes of the W boson exchange graphs with their explicit calculation 
of counterterms in the MS scheme a logarithm ln/x 2 /m 2 remains left. 

One can observe that for the QCD corrections the coefficients of the In /x 2 /m 2 
and the lnM^/m 2 terms in eq.(^) are identically the same as in the pure elec- 
troweak result of eq.flj^). Thus these logarithms are characterized by the same 
correction factor (1 + c^/V) as it is known from the pure QCD corrections to the 
Born decay rate. With the electroweak result being well established in the litera- 
ture, the explicit /i dependence in the OS result eq.® is known to be compensated 
by the running electroweak running coupling constant, thus giving a RG invari- 
ant result. The common correction factor (l + a s /7r) therefore also guaranties the 
/^-cancellation in our 0(aa s \nm^) result. 

Having gained more experience and insight in the heavy mass expansion as a 
practical approach of an effective theory since our previous work |7|, let us add 
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an additional comment concerning the appropriate choice of scale as argument 
for the running mass and coupling constant. 

The leading 0(aa s m 2 ) calculation exhibits a remarkable structure. All in- 
tegrals factorize and can be classified into two cathegories, both of which are 
separately finite and gauge invariant. The first one is characterized by a factor- 
ization into a two loop massive tadpole integral and a massless one loop p- integral, 
where the gluon is part of the former. Since this integration is affected by only 
one mass scale namely the top quark mass, it is instructive to evaluate a s at the 
scale mf for these contributions. Contrary to this the second cathegory comprises 
all corrections which factorize into a one loop massive tadpole and a two loop 
massless p-integral, with the gluon lines contained in the latter. The only scale 
in the massless p-integral is the energy regime of the process under consideration 
and it seems to be appropriate to calculate a s at this scale. 

The result for the 0(aa s m 2 ) term can now be written as a factorized form, 
where both scales are separated and higher order QCD corrections of order a 2 
are neglected. 

Recalling that pure QCD corrections to T(Z — > bb) are given by (1 + a s (s)/n) 
the following interpretation is at hand. The virtual top quark is of purely elec- 
troweak origin and its effect can be accounted for by effective vertices. QCD 
corrections enter twofold. On the one hand they may be part of the effective ver- 
tex with the relevant scale being m 2 . On the other hand these QCD corrected 
vertices are dressed by a virtual gluon, thus being multiplied by the correction 
factor (1 + a s (s)/n). The scale of the top mass is in both cases p 2 = m 2 because 
the overall mf factor always results from the tadpole integral. We have not found 
such a factorization in the \nm 2 term, mainly because the two classes are not 
separately finite. 

To summrize all contributions discussed in this paper, we use the common 
parameters p = 1 + 5p and k — 1 + 5k through which the Z decay rate into two 
b quarks can be expressed: 



T(Z -> bb) 



G F M\ 




w 




6ft 



non— univ. 



Gp 



2y/2n 2 




(12) 



+^Mf In 



7T 
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w 



) 



6 



+ Ar 

r^_non— univ. . _ jj^r non— univ. 

2 P 



OS 

rem 



Figure 2: Leading order and sum of leading and next-to-leading order QCD 
corrections. 

The leading mf QCD correction is positive in the 05* and negative in the MS 
scheme. In both cases the corresponding lnmf term has the opposite sign. The 
leading order OS result is reduced by about a factor one half. The modulus of 
the next-to-leading term in the MS scheme is slightly bigger than the term, 
so that the sum is small and positive. 

In Figure the leading order and the sum of the leading and next-to-leading 
order corrections are plotted against the top mass. It is easy to see that the 
inclusion of the next-to-leading correction reduces the difference of the predictions 
between both schemes considerably. 

Figure |3] compares the QCD corrections with the pure electroweak result both 
in the MS and in the OS scheme. One can recognize that in the MS scheme the 
corrections are numerically tiny (2.5 • 10 -4 ). 



Figure 3: Pure electroweak corrections (solid line), electroweak + QCD correc- 
tions (dashed line). The upper curves belong to the MS scheme. 



To conclude, we calculated QCD corrections to the known electroweak result 
for the partial width Y[Z — > bb) in the limit of a heavy top quark. We used 
an expansion in the inverse top mass and calculated the next-to-leading term of 
order 0(aa s lnmf). 
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